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1. Introduction 

In a recent article [IJ (hereafter referred to as I) the finite-size scaling behaviour 
of a three-dimensional system with a film geometry L x oo 2 has been investigated 
within the mean spherical model with Neumann- Neumann and Neumann-Dirichlet 
boundary conditions and surface fields hi and h L acting at the boundaries. The 
obtained results imply the well known exponent A° = 1/2 for the ordinary surface 
phase transition at a Dirichlet boundary, and the emergence of a new critical exponent 
Af 3 = 3/2, characterizing the Neumann boundary (for a general review on surface 
critical phenomena see, e.g. [0-0], and on finite-size scaling 0, || - 0). The 
conjecture has been made that the latter critical exponent corresponds to the special 
(surface-bulk) phase transition within the model. The last is in consistence with the 
general expectation for the finite-size scaling form of the free energy for this type of 
phase transitions if one assumes that the crossover exponent $ = 0, as it is for three- 
dimensional 0(n) models JJ. It has also been derived that the critical exponent of 
the local surface susceptibility xi,i is 7i b i — 1- The same result is known to hold for 
the spherical model with enhanced surface couplings under Dirichlet-Dirichlet boundary 
conditions 0. Unfortunately, in that case the model quite unphysically predicts that 
the surface orders for a sufficiently large enhancement at some temperature above the 
bulk critical one even for d = 3. This is no more the case when one improves the 
model by introducing a second spherical constraint on the spins at the boundaries ||, 
since the only critical point that remains for d < 3 is the bulk one. Then for d — 3 
the exponent j° t = —1 corresponds to an ordinary phase transition || 10]. In I the 



case of equal bulk and surface couplings was considered and the question if and how 
the surface behaviour of the system with Neumann-Neumann boundary conditions will 
change under additional spherical constraints on the spins at and near the surfaces was 
left open. One of the aims of the present article is to contribute in clarifying that point. 
To this end we consider the critical behaviour of the local surface susceptibility 



Xi,i(T; p) = Km -Ld 2 f L (T, hi, h L ; p)/dh{ 



(1.1) 

\h 1= h L =0 v ' 



in the case when the mean square value of the spins at the boundaries is fixed at some 
positive number p by an additional spherical constraint. The model defined in this 



way will be termed modified spherical model. In equation (|1 . 1|) and in the remainder 
/l(T, hi, hL', p) denotes the free energy density (per ksT and per spin) of a three- 
dimensional hypercubic lattice system with a film geometry L x oo 2 at temperature 
T. Neumann-Neumann boundary conditions are imposed across the finite dimension of 
extent L. Surface fields hi and h L are supposed to act at the surfaces bounding the 
system. Since we are interested in the case hi = hL, only one additional constraint 
on the boundary spins is imposed. It turns out that the behaviour of i n the 
vicinity of the bulk critical point T c depends crucially on p. It will be shown that 
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only if p = p c := 2 — (12i^ c ) _1 = 1.34053. . . one obtains the previously found value 
7i ; i = 7* = 1. If p < p c then xi,i has the singularity characteristic of the spherical 
model with Dirichlet boundary conditions, i.e. 71,1 = 7° 1 = —1. When p > p c then Xi,i 
diverges exponentially as T — > T c + , which reminds the behaviour of a two-dimensional 
0(n), n > 2, model close to T = 0. The calculation of the mean square value p s of the 
spins at a Neumann boundary in the standard spherical model elucidates the appearance 
of the critical value of p = p c : it turns out that p c = p s . Moreover, we note that the 
second spherical field, to be denoted by v (see equation ( |2.1| ) below), can be considered 
as a free parameter. It will be shown that changing it one interpolates continuously from 
Neumann, via mixed, to Dirichlet boundary conditions. Following I, under Neumann 
boundary conditions we mean here the case when the interaction of the finite system 
with the "environment" is modelled by letting the spins surrounding the system take 
the same values as their nearest neighbour inside the system. Under Dirichlet boundary 
conditions this interaction is modelled by fixing to zero value the spin configuration 
outside the system. (For a precise mathematical definition of the boundary conditions 
see I.) The mixed boundary conditions correspond then to the situation when the spins 
surrounding the system are set to take values proportional (but not equal) to those of 
their nearest neighbour inside the system. Obviously, the above terminology is justified 
by analogy with the continuum limit. Note that for any v , just due to the symmetry 
which arises from the identical boundary conditions and fields (hi = hi) at the opposite 
surfaces, the system models by itself an analogue of a Neumann boundary at the middle 
layers. Therefore, if one considers the local surface succeptibility xi,i f° r the /-th layer, 
one would expect to obtain the critical exponent for the Neumann boundary, 7^ = 1, 
for I around the middle of the system. The last is obviously true even if the system 
is with otherwise Dirichlet boundary conditions. Finally, it will be shown that if v is 
a given function of the temperature, one obtains 7* = 2, which is the corresponding 
value for the 0(n) model in the limit n — > 00. 

As it is well known, the infinite translational invariant spherical model is equivalent 
to the n — > 00 limit of a similar system of n-component vectors [[□]] , |12] . However, 
the spherical model with surfaces (or, more generally, without translation-invariant 
symmetry) is in fact not such a limit |13] (for the results available for the spherical 
model see, e.g.,@, ||, @, an d references therein). In other words, the spherical model 
under nonperiodic boundary conditions is not in the same surface universality class as 
the corresponding 0(n) model in the limit n — > 00, in contrast with the bulk universality 
classes. The last becomes apparent when one investigates surface phase transitions for 
an 0(n) model in the limit n — > 00. In that case one obtains H A x = l/(d — 2) 
(i.e. Ai = 1 for d = 3) for ordinary and Ai = 2/(d - 2) (i.e. A x = 2 for d = 3) 
for special phase transitions. It is believed that the corresponding equivalence will be 
recovered if one imposes spherical constraints in a way which ensures that the mean 
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square value of each spin of the system is the same |13j (unfortunately such a model is 
rather untractable). One of the goals of the present work is to see if, and up to what 
extent, the behaviour of the system with two spherical fields will be closer to the 0(n) 
model in the limit n — > oo, in comparison with the standard mean spherical model (with 
only one spherical field). 

The article is organized as follows. In Section 2 we describe the model and present 
convenient starting expressions for the mean spherical constraints, the free energy 
density and the local surface susceptibility. Our main results on the behaviour of Xi,i 
and xi,i as a function on T and p are given in Section 3. The paper closes with a short 
discussion given in Section 4. 



2. The model 



We consider explicitly the three dimensional mean spherical model with nearest- 
neighbour ferromagnetic interactions on a simple cubic lattice. At each lattice site 
r = (^i,r 2 ,r 3 ) G Z 3 there is a random (spin) variable a(r) G R and the energy of a 
configuration <j a = {a(f),f G A} in a finite domain A C Z 3 , A = L\ x L 2 X L3, containing 
|A| sites, is given by 

W^KIK, h A ; s) = -Ka\ ■ ■ a A + sa\ ■ a A + wj • a. - h\ ■ a A . (2.1) 
Here the |A| x |A| interaction matrix Q^' can be written as 

= (AS ri) + 2EJ x (A { P ] + 2E 2 ) x (A^ + 2E 3 ), (2.2) 

(t) 

where x denotes the outer product of the corresponding matrices, A- is the Li x Li 
discrete Laplacian under boundary conditions Tj, and Ei is the Lj x Lj unit matrix. 
In equation ( [2.1| ) j3 = 1/ksT is the inverse temperature; K = j3J is the dimensionless 
coupling constant; h A = {h(r),f 6 A}, with h(r) G R, is an external magnetic field; 
s and v are the spherical fields which are to be determined from the mean spherical 
constraints (see below); a s = {a(f),f G S}, S = {(r±, r 2 , l)U(ri, r 2 , L 3 )}, r% — 1, . . . , Li, 
r 2 = 1, . . . ,L 2 . 

The free-energy density of the modified mean shperical model in a finite region A is 
given by the Legendre transformation 

PfW(K, h A -p) := sup {-\K\- 1 In Z£\K, h A ; s,v)-s- pv\S\/\A\\ , (2.3) 

s,v 1 J 

where \S\ is the total number of spins at the boundaries S and 

Z^{K,h A ;s,v) = J exp [-^\a A \K, h A ; s, v)\ ]J da(r) (2.4) 

is the partition function. The supremum is attained at the solutions of the mean 
spherical constraints 

(4 ■ O = |A| (2.5) 
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and 

(at ■ a s ) = p\S\, (2.6) 

where (• • ■) denotes expectation value with respect to the hamiltonian j3TC^'(a A \K, h A ;s). 
Let us denote by —2 + 2cos<£>2 .(ki), h — 1, . . . , L i: i — 1,2,3, the eigenvalues of the 
matrix A\ T ^ . Let us further suppose, say, periodic boundary conditions across L\ and 
L 2 and Neumann- Neumann boundary conditions across L 3 . Then, by direct evaluation 
of the integrals in the partition function {2A), after taking the limit Li,L 2 — > oo at a 
fixed L3 = L, one obtains for the free energy 



1 K 

-log--6K + 

l 7T 



su p(t^ ^2 Elog[0 + 2^(1- cos ^)] + 

<j>, w { L 87T 2 Jo Jo ^[ ^ 



2[l-cosy>2(Jfc;w)] 



AKL 



E 

k=l 



1 


i=l 




h£\k;u)\ 2 


+ 2 


1 — cos<^2(fc; to)} 



K 
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Here —2 + 2 cos <^2(^> cu), = 1, . . . , L are the eigenvalues of the matrix 



A 



(n), 



A[ n) - u (6 ltl + 6 LlL ), 



and 



(2.7) 

(2.8) 
(2.9) 



h^\k]uo) = h lU n L (l,k;uj) + h L u n L (L,k;Lu) 

where {w"(r, k; u>), r = 1,...,L}, k = 1,...,L are its eigenvectors; the superscript n 
stays there for Neumann-Neumann boundary conditions. In equations (2.7) - fl2.9p use 
has been made of the following definitions 

<f) = s/K-6, (2.10) 

and 

u = v/K. (2.11) 
From the requirement for existence of the partition function one has the constraint 

+ 2 min [1 - cos <pl(k; cu)} > 0. (2.12) 

k=l,...,L 

The eigenvalues and the eigenvectors of the matrix A^ (u;) can be obtained in a way 
simmilar to the one used in || and [[14]]. The results are: 

i) For given L and u, when |1 — uo\ 7^ 1, the numbers (p r l(k;uj),k = l,...,L are the 
L roots of the equations 

1 - ^ = ^ — i (2.13) 



sin 






sin 




1)<P 



6 



and 



1-iO 



cos 




COS 





(2.14) 



with <Re(y?) < tt and Im(<£>) > 0. For concreteness and simplification of the notations 
below, without loss of generality in the final results, we assume L to be an odd integer. 
Then, it is easy to see that equation ( |2.13| ) possesses (L— l)/2 solutions of the specified 
type, whereas equation ( |2.14j ) gives the remaining (L + l)/2 solutions. From ( |2.13| ) and 
( |2.14| ) one obtains that (for fixed L and k ) 

= I 2sin 1 (2 15) 

du L(1-cj) 2 -2(1-cj)cos^ + 1 + L- 1 [1-(1-u;) 2 ]' V ' ; 

Further, if |l — u\ < l all the L roots are real. In that case there is only one root of 
equation flTl3| ) per interval (2^k/{L - 1), 2n(k + 1)/(L - 1)), k = 0, . . . , (L - 3)/2. 
Similarly, equation Q2.14j ) has only one root per interval (it(2k — 1)/(L — 1), 7r(2fc + 
1)/(L — 1)), k — 1, . . . , (L — 3)/2, and one root in each of the intervals (0, tt/(L — 1)) 
and (tt — tt/(L — 1), it). Let us consider now the case |1 — u\ > 1. Then, if uj < 0, one 
again has only one root of equation ( |2.13| ) per interval [2i\kj (L — 1), 2n(k + 1)/(L — 1)), 
k = 1, . . . , (L — 3)/2, and, similarly, one root of equation ( |2.14j) per interval (ir(2k — 
l)/(L - 1), n(2k + 1)/(L - 1)), k = 1, . . . , (L - 3)/2, and one root in (tt - ir/(L — l),7r), 
i.e. altogether L — 2 real roots in the interval (0, tt). The remaining two roots are given 
by 

<p = i log(l - uj) ± O ((1 - . (2.16) 

(Strictly speaking the roots are degenerate only up to exponentially small corrections.) 
In the case to > 2 one again has L — 2 real roots in the interval (0, tt) and the remaining 
two roots are then given by y?o = tt + i log(u — 1) ± O (^(u — 1)~( L-1 )) . 

ii) The components of the eigenvectors {w™(r, k;u),r = 1, . . . , L}, k = 1, . . . , L of 
the matrix A^ 1 (a;) are given by the expression (|1 — u\ ^ 1) 

n (rh a [2 sin[ry£(fc; ^)] - (1 - ^) sin[(r - cj)] 

ik ' ' J V £ {(1 - uj) 2 - 2(1 - u) coscpt(k;co) + 1 + L^[l - (1 - ^j) 2 ]} 1 ^ 

iii) For completenes we give also the results for the well known case |1 — u\ = 1 
(see, e.g, I, 0). Then cpl(k;0) = vr(A; - 1)/L, yg( fc; 2) = irk/ L, k = 1, . . . , L, and 
the components of the eigenvectors are u n (r, A;; 0) = J (2 — 5fe i i)/Lcos[(r — l/2)(f 1 l(k; 0)] 



and u™(r,k;2) = J2/Lsin[(r — l/2)<p2(k; 2)], respectively. 

Finally, we remind that we are interested mainly in the behaviour of the local surface 
susceptibility for which we obtain from equations ( |1 . 1|) and (2.7) 

* l( T; p) = -L Mm V M' 1 -*""' (2.18) 
A1,u 2iiri-oo^0 + 2[l - cos <fl(k;cu)} K ' 
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If, instead of the local susceptibility at the surface of the system, one is interested in 
the local susceptibility of the l-th layer, Xi,i(T; p), the corresponding result reads 

Xi,i(T; P) = 7^7 I™ E J,i 7i 2 - 19 

2a l-+oo J + 2[1 — cos <^2(k; u;)J 

The above expression can be obtained in a way analogous to the derivation of xi,i by 
imposing a local magnetic field hi on the spins in the l-th layer. 

To determine the behaviour of the spherical fields and to one has to analyze 
equations fl2J| and (glU) . From @, (2.7), ( CT) , (PH) and ([2151) one obtains 
explicitly the set of equations 



2K = — d6 x d9 2 -Y, { <j> + 2^{1- cos ei) + 2[l- cos cpl{k;u))\ (2.20) 
4tH Jo Jo L k=1 { . =1 J 

and 

2Kp = -L /** d0i /^ d0 2 y£2 sin 2 ^(fc; <*>) x 

47T Z JO JO L 

L) + 2J2(1 -cos6i) + 2[1 -cos^2(fc;^)]| x 
(1 - cj) 2 - 2(1 - w) cos(fl(k; oj) + 1 + L -1 [l - (1 - a/) 2 ]}" 1 . (2.21) 



These equations determine the point at which the finite-size free energy density (2.7), 
which is analytical and strictly concave function of and to in the domain given by 
inequality ( |2.12|) , reaches its global maximum. Clearly, in the thermodynamic limit the 
free energy density is independant of the surface spherical field to. As it is well known, 
for all K > K c its supremum sticks to the endpoint 0o — of the allowed interval 
0o > 0, where the bulk free energy density is finite. When K < K c , the supremum is 
attained at a point O = <j>o(K) > 0, which satisfies the limit form of equation (|2.20| ) 

m, 



2K = W 3 ((P ), (2.22) 



where 



f"K 


P7T 


\ dd x -- 


■ / dd d 


lo 


Jo 



+ 2J2(1 -cosOi 



(2.23) 



is the (i-dimensional Watson integral, Wd(0) = 2K C . In general, the solutions and to of 
equations ( |2.20| ) and (2.21) can be written in the form = O + A0, to = io + Ato, where 
A0 and Ato tend to zero when L — > oo, and 0o and too are solutions of the corresponding 
equations where the limit L — > oo is taken. 

Equations ( |2.16| ) - ( |2.23| ) provide the basis for our further analysis. Before passing 
to it we note that, instead of considering to as a variable that has to be determined 
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from equations ( |2.20| ) and (2.21), one can consider it as an additional free parameter. 
Then, from equation (|2.8|) it is clear that u = yields the standard spherical model 
with Neumann-Neumann boundary conditions, whereas to = 1 yields the same model 
under Dirichlet-Dirichlet boundary conditions. When < uj < 1 we have mixed (or 
"intermediate" |To| ) boundary conditions which interpolate between the above two 
extreme cases. Therefore, in this way one should reproduce the previously known 
results for the properties of the local susceptibilities. In addition, as we shall see later, 
by choosing u to be a given function of the temperature, one can define an effective 
spherical model with 7^1 = 2, which corresponds to the critical exponent for the surface- 
bulk phase transition within the 0(n) model in the limit n —>■ 00. 



3. Critical behaviour of the local susceptibilities 



We study here the critical behaviour of the local susceptibilities Xi,i an d Xi,i f° r I close 
to the middle of the system. 

From equations ( |2.18|) and ( |2.17| ) we obtain for the surface susceptibility 



xi,i( t ;p) 



— lim — 



sin 2 <£>£(&; u) 



+ 2[1 - cos (fl(k;uj)] 
1 - 2(1 - uj) cosip n L (k; u) + (l- uj) 2 + [1 
where the limit L 



{l-ujf]/L 



(3.1) 



00 in (3.1) is to be taken over the finite-size solutions uj and <p of 
equations (|2.20|) and (2.21). If |1 — ui\ < 1, from the properties of <^£(fc; uj), k = 1, . . . , L, 
described in Section 2, it follows that as L — > 00 the sum in equation (3.1) tends to the 
corresponding well defined integral 

K n Jo 



xi,i( t ;p) 



sin 2 ip 



+ 2(1 — COS<y9)] [1 — 2(1 — U!q) COS if + (1 — LOq) 



(3.2) 



The integral can be taken exactly [16] with the result 

xi,i( t ;p) 



K 



yJ</>0 (4 + O ) + 00 + 2cj 



(3.3) 



When |1 — ui\ > 1 one has to take into account the contribution of the two complex 
roots which turns out to be of the same order as the contribution of all other roots. 
The contribution of the latter L — 2 roots is again given by the integral in the r.h.s 
of equation (|3.2|) . Performing the calculations one ends up with the same analytical 
expression for Xi.iC^i P) as ^ ne one gi yen by equation (|3~3"1). 

The surface spherical field u Q satisfies the corresponding limit form of the spherical 
constraint (2.21) at fixed 0o — for K > K c , and 0o — 4>o{K) for K < K c . The 
right-hand side of this equation can be treated in a way similar to that for (3.1). When 
|1 — uj\ < 1, due to the properties of the roots ip^k^u)), k = 1, . . . L, the sum in (2.21) 
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converges as L — > oo to the corresponding well defined integral, which can be taken 
exactly. Performing this procedure, one obtains finally 



2Kp = G 3 (<p ,u ), 



(3.4) 



where 



PIT 




/ d0!.- 
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d-l 



d-l 



,1/2 



+ 2^(1 - cos 6i 



i=l 



+ 4 + 2^(l-cos^ 



i=l 



1/2- 



(3.5) 



When 1 1 — oj \ > 1 , one has to treat separately the contribution from the two complex 
roots. The contribution from the L — 2 real roots leads again to a well defined integral 
that can be taken exactly. As for xi.i{T;p), the final result is given by the same 
analytical expression as in the case |1 — uj\ < 1, i.e. equation (|3.4j) is actually valid 



for all uq (the restrictions on ujq and 0o stemming from the constraint (|2.12|) are stated 
below) . 

Let us denote by G 3 ((f>,uj) the branch of the function G 3 (<p,u) defined for uj > 
and by G 3 (<p,u) the one for uj < 0. Then, by means of identical transformations it is 
easy to show that 



G 3 ((j),uj) 



2 n + 
3 



\0J\ 



uj{2-uj) 



Wo 



uj 



l-UJ 



and 

Gj(<j ) ,u) = (l-u)- 1 



2WM- \ + \<t>wM 

o o 



(3.6) 



(3.7) 



UJ 



[1 -uj)(2-uj)ti 2 Jo 



de 1 



uj 



2-uj 



+ 2^(1 -cosflj 



&e 2 

.1/2 



2^(1 - cos e t 



i=l 



i=l 



+ 4 + 2^(1 -cos^ 



.1/2- 



i=l 



Finally, in the limit L — > oo the constraint Q2.12|) for the existence of the partition 
function yields the allowed domain of values of the spherical fields, 



> 



(3.f 



0, if uj > 

uj\I (1 - u ), if uj < 0. 

From equation (|3.3| ) it follows that the above inequalities imply xi,i(T'i p) > 0, as it 
should be expected on general physical grounds. 

As it is evident from equation (3.5), G 3 ((p, uj) is a monotonically decreasing function 
of uj which tends to zero from above as uj — * +oo. Due to inequalities (|3.8|), at <fi = we 
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have to consider it on the half-line to > 0, where it is bounded from above by its value 
at to = 0, see equation (3.7), 

G 3 (0, 0) = 2K C - 1/6 := 2K cPc . (3.9) 

On the other hand, if > 0, the definition domain of G^{(j),Lo) is restricted by ( p.8|) to 
the half-line to > u>i(4>), where 

^ 1 (0) = -(0 + 2 /4) 1 / 2 -0/2. (3.10) 

From the representation ( |3.6| ) and the known expansion of Wzipc) as x j 0, 

W 2 (x) = (47r)- 1 lnx- 1 + 0(l), (3.11) 

it follows that G 3 {4>,lo) diverges logarithmically to +oo as to | u>i((f>). 

Before passing to the analysis of the above equations, in order to determine the 
behaviour of xi,i(T;p), let us first consider the simpler case of to as a free parameter. 
Then, for Neumann- Neumann boundary conditions one has (see equations (|2.2| ) and 
Q2.8Q ) to — 0, whereas to — 1 for Dirichlet-Dirichlet boundary conditions. Thus, from 
( p.3|) and the well known behavior of 0o m the vicinity of the bulk critical temperature 
0o — [8ir(K c — K)] 2 [0, one immediately obtains all previously known results for the 
critical behaviour of the local surface susceptibility 0, [|Hj: 

a) Neumann-Neumann boundary conditions {to = 0; the result given below follows 
directly from equation (3.5) in for h\ = h£) 

Xi,i(T) = (2K)- 1 {0 o /2 + [0 O (1 + o /4)] 1/2 p , (3.12) 
i.e. 71,1 = 1, and 

b) Dirichlet-Dirichlet boundary conditions (u = 1; see equation (61) in [ [jl| ) 

Xi,i(T) = (2K)' 1 {l + o /2 + [0o (1 + 0o/4)] 1/2 }^ 1 , (3.13) 

i.e. 71,1 = -1. 

For to 7^ 0, 1 one has the case of the so-called intermediate ]TS[ boundary conditions. 
As it is clear from ( |3.3| ), xi,i diverges in the vicinity of T = T c if and only if to = 0, i.e. 
under Neumann-Neumann boundary conditions. 

Let us now comment on the critical value p c of the parameter p, defined in equation 
Q3.9| ). By using translation invariance argument, for the mean square length of the spins 
at the Newmann boundary of the standard spherical model with one global spherical 
field one obtains in zero magnetic field 

£ , iffiM --, ™ (3.14) 



fci,fc 2 ,fc3 = l 



+ 2ELJ1 - cos(2 7 rA; l /L i )] + 2[1 - cos(7r(fc 3 - 1)/L 3 )] ' 
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In the limit of an infinite film geometry this equation yields (L 3 = L is kept finite) 

lim (<7 2 (ri,r 2) l)> = hw z {(j>) - 1/12 + 0W 3 (</»)/12 + W 2 (0)/2L]. (3.15) 

Hence, at the critical point if = A" c of the infinite system (L = oo), by taking into 
account the bulk spherical constraint at 0o = 0, namely H^O) = 2K C , one obtains 
that the mean square length p s of the spins at the Neumann boundary of the standard 
shperical model equals precisely the critical value p c for the surface spins in the modified 
spherical model. 

Finally, we note that by taking 

u = -8n(K c -K), (3.16) 

one obtains for the considered system with layer geometry and Neumann-Neumann 
boundary conditions 7^1 = 2, which is the corresponding critical exponent for the 0(n) 
models in the limit n — > 00. In that case 7J 1 exists too, and 7^ 1 — 1. Obviously, such 
a choice of u defines an effective Hamiltonian that leads to an exactly solvable model 
with the critical exponents stated above. 

Now we pass to the analysis of the behaviour of Xi,i(T; p) given by equation ( |3.3| ) 
where uq is determined as a function of K and p from equation (|3.4| ). 



3.1. Critical behavior of the local surface susceptibility 

Here we confine our analysis to the surface critical regimes that emerge on approaching 
the bulk critical temperature from above, i.e. when K = K c + AK, where AK < and 
I AK\ — > 0. Then, as it is well known, the leading asymptotic form of the bulk spherical 
field follows from the asymptotic expansion 

H/ 3 (0) = 2K C - (47r)-y /2 + 0(0), I (3.17) 

and reads JjlJ 

= QAtt 2 \AK\ 2 , AK t 0. (3.18) 

From expression ( |3.3|) it is clear that the local surface susceptibility may exhibit 
divergent behaviour in two different regimes: (a) when w | 0, and (b) when uq j 
uj\{4>) t 0. As it is clear from equation ( |3.4j ) and the above mentioned properties of 
the function G 3 (0, cu), the first regime may occur only when Kp j K c p c , which, in view 
of our assumption AK j 0, requires p = p c . The second divergent regime of the local 
surface susceptibility takes place at any fixed p > p c . Below we derive the leading-order 
asymptotic solutions for co> in each of the two cases. 

Case (a): p = p c . To obtain an asymptotic expansion of Gs(4>, u) in both arguments 
I and uj I 0, we notice that the integral in the r.h.s of equation (3.7) diverges at 
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<p = <jj = and the divergence arises from the integration over the neighbourhood of the 
point 61 = 62 = 0. Therefore, its leading-order asymptotic behaviour is given by the 
small argument expansion of the trigonometric functions which yields 

G+(0, u) = 2K cPc - (27r)- V /2 + (w/2tt) ln(0 1/2 + u) + 0(0) + 0(u). (3.19) 

By setting p = p c , we obtain that co>o I obeys the asymptotic equation (AK | 0) 

-(w q /2tc) ln(87r|AK| + u ) = \AK\/(6K C ). (3.20) 

The solution which tends to zero from above as \AK\ — > is 



7T 


AK 




3X c ln(87r| 


AK\) 



^0 - - or ^ : ' ,' (3.21) 



Obviously, this critical regime leads to 71,1 = 1. 

Case (b): p > p c . The asymptotic behaviour of Gs((j), uj) as <fi | and | 0, so that 
uj > is readily obtained from the exact representation (|3.6|) and the expansions 

( grp and (|3T9|) : 

G?a (0, w) = 2K c p c - (2tt)-V 1/2 - (|w|/2tt) ln^ 1 / 2 - + 0(0) + O(cj). (3.22) 
At fixed Ap > the leading-order equation for the surface spherical field becomes 

-(|wo|/27r) ln(87r| AK\ - \u \) = 2K c Ap. (3.23) 
Assuming |o; | = 87r|Ai^| — x, where x = o(\AK\), one obtains 

uj ~ -8tt| AK\ + exp (-^^"J ■ (3-24) 

Therefore, in this critical regime the local surface susceptibility diverges exponentially 
as the bulk critical temperature is approached from above: 

XM(T;p)^exp(^). (3.25) 

This behaviour reminds the one of a two dimensional system close to T = 0. The fact 
that the surface is coupled to an infinite three dimensional system is reflected in the 
replacement of T = by the bulk critical temperature T = T c . 

Finally, if p < p c it is easy to see that (|3.4| ) has a finite solution u (p,K), where 
< Uq < 1, when AK f 0. The last actually follows from the inequalities 

G 3 (O,0)>^3(0)>O 3 (1,0). (3.26) 

Thus, if p = 1 the local surface susceptibility Xi,i(T c ; 1) is finite. 



13 



3.2. Critical behavior of the local susceptibility around the middle of the system 

For the local susceptibility Xi,i(T; p) from ([2.171) an d Q2.19|) one obtain explicitly 

1 1 * {sinK(fc; u)] - (1 - u) sin[(l - l)yg(fc; a;)]} 2 
X M (i , pj - — lim - ^ ; rr: — - X 



Kl-*ooL^ x (j) + 2[1- cos (fl(k;u)] 

[l - 2(1 - to) cos<pl(k; u) + (l- uf + [1 - (1 - u) 2 ]/ ' l)' 1 . (3.27) 



We will be interested only in the behavior of this quantity around the middle of the 
system. Let us set I = (L + l)/2. Then, from (3.27) it follows 

Xi,i(T; p) = ^ lim j ]T [0 + 2[1 - cos <p n L {k; iv)]]' 1 x 

\ L [1 - 2(1 - u) costpKk; ui) + (1 - uof] + 1 - (1 - uof J ' 1 ; 

where the summation is over the roots of equation (|2.14j) only. Having in mind the 
properties of the roots </?2(fc; ui), it is easy to see that in the limit L — > oo this equation 
leads to 

Xoo,oo(T;p) = ^^ 1 (0o) (3.29) 

for any u. We recall now that considering u as a free parameter, at uo = one has 
the standard spherical model under Neumann-Neumann boundary condition and at 
uj = 1 the corresponding one with Dirichlet-Dirichlet boundary conditions. The above 
result shows that the behavior of Xoo,oo(^; p) does not actually depend on the boundary 
conditions. From the temperature dependence of 0o around the bulk critical temperature 
0o ~ [8tt(K c — K)] 2 ]T7[ and the expansion of Wi(<f)) for small values of the argument 

W 1 {<P) = \<P- 1I2 + 0{<P 1 '% (3.30) 

we obtain 7oo,oo = 1- It is clear that the same will be true for any layer at finite distance 
from the middle of the system. The above result has been derived in |TIJ for a spherical 
model under Dirichlet-Dirichlet boundary conditions (see Eq. (82) in ||10|| ). Here we 
simply show that it does not depend on the boundary conditions if they are identical at 
both the boundaries: just due to the symmetry the system models by itself an analogue 
of the Neumann boundary at the middle layers. 



4. Discussion 



In the present article the surface critical behaviour of a modified three dimensional 
mean spherical model with a L-layer film geometry under Neumann-Neumann boundary 
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conditions is considered. The standard spherical model is modified in the sense that 
in addition to the usual bulk spherical constraint a second spherical field is included in 
the Hamiltonian to fix the mean square value of the spins at the boundaries at some 
value p > 0. We are interested mainly in the upper critical behaviour of the local 
susceptibilities an d Xl,i with I close to the middle of the system. The surface 
susceptibility xi,i an d the local susceptibility Xoo,oo are evaluated exactly and the 
corresponding results are given by equations ( |3.3D and ( |3.29| ), respectively. 

It is shown that the behaviour of Xi,i(T; p) depends crucially on p. At p = 1 we find 
that xi,i is finite at the bulk critical temperature T c , in contrast with the recently derived 
value 7i i = 1 in the case of just one global spherical constraint. The result 7^1 = 1 
is recovered only if p = p c = 2 — (12i^ c ) _1 , where K c is the dimensionless critical 
coupling. When p > p c , the local surface susceptibility xi,i diverges exponentially as 
T — > T c + , see equation ( |3.25| ). The calculation of the mean square value p s of the spins 
at the Neumann boundary in the standard spherical model elucidates the appearance 
of the critical value of p = p c : it turns out that at the bulk critical point p c = p s , see 
equation Q3.15Q at K — K c , = and L = 00. As it is expected, the behaviour of 
the local susceptibility Xoo.oo turns out to be independent of the boundary conditions if 
they are the same at both boundaries. Just due to the symmetry, the system models 
by itself an analogue of the Neumann boundary at the middle layers which leads to 
7oo,oo = 1 (see Section 3.2 for details). By considereing the second spherical field as an 
independent free parameter, we rederive in an uniform way the previously known critical 
properties of the local surface susceptibility. They follow directly from equation ( |3.3j ) at 
uj = 0, for Neumann- Neumann (see ( |3.12| )), and u — 1, for Dirichlet-Dirichlet boundary 
conditions. For uj 7^ 0, 1 equation ( |3.3| ) gives the corresponding result for the so-called 
"intermediate" Jl5| boundary conditions. From these results one concludes that 
diverges at T = T c only under Neumann boundary conditions. Finally, an effective 
hamiltonian which leads to an exactly solvable model with 7^ = 2, the value for the 
n —>■ 00 limit of the corresponding 0(n) model, is proposed. It is given by equation 
( |2.1|) where one has to set v = —8iiK(K c — K), see ( 13.161) . 

We emphasize that the spherical model under nonperiodic boundary conditions is 
not in the same surface universality class as the corresponding 0(n) model in the limit 
n — > 00, in contrast with the bulk universality classes. For example A° = 1 and Af 5 = 2 
for the 0(oo) model, but A° = 1/2 and A^ b = 3/2 for the spherical model. The results 
presented above show that the properties of the model are improved by introducing a 
second spherical constraint in the sense that they are closer, in a certain way, to the 
corresponding ones for the 0(n), n — > 00, model. It seems clear that in order to obtain 
"correct" surface critical properties, one has to impose a separate spherical constraint 
on each layer paralell to the surface. 
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